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Abstract 

We propose two new simple lepton flavor models in the framework of the S '4 flavor symme¬ 
try. The nentrino mass matrices, which are given by two complex parameters, lead to the 
inverted mass hierarchy. The charged lepton mass matrix has the 1-2 lepton flavor mixing, 
which gives the non-vanishing reactor angle @ 13 . These models predict the Dirac phase 
and the Majorana phases, which are testable in the future experiments. The predicted 
magnitudes of the effective neutrino mass for the neutrino-less double beta decay are in 
the regions as 32 meV < \mee\ ^ 49 meV and 34 meV < \mee\ ^ 59 meV, respectively. 
These values are close to the expected reaches of the coming experiments. The total sum 
of the neutrino masses are predicted in both models as 0.0952 eV < < 0.101 eV 

and 0.150 eV < < 0.160 eV, respectively. 


*E-mail address: shimizu@muse.sc.niigata-u.ac.jp 
lE-mail address: tanimoto@muse.sc.niigata-u.ac.jp 



1 Introduction 


The neutrino experiments are going on the new step to reveal the CP violations in the 
lepton sector. If the neutrinos are the Majorana particles, there are one Dirac phase and 
two Majorana phases PEIE], which are sources of the CP violation. The T2K experiment, 
which has already confirmed the neutrino oscillation in the —)• z/g appearance events [1], 
provides us the new information of the CP violating Dirac phase by combining the data of 
reactor experiments 0-0. On the other hand, the neutrino-less double beta decay {0h'(3(3) 
experiments lower gradually the upper-bound of the effective neutrino mass ruee around 
(P(IOO) meV. If the neutrino mass hierarchy is the inverted one, one can expect to observe 
the Oz//3/3 in the near future. Its magnitude depends on the CP violating Majorana phases. 
Thus, the three CP violating phases can be observed in progressive neutrino experiments. 
Those CP violating phases provide us the important information to test the flavor model of 
neutrinos. There appear many works to test models with the CP violating Dirac and the 
Majorana phases of neutrinos in addition to the mixing angle sum rules [9]-|32]. 

The recent experimental data of neutrinos give us a big hint of the flavor symmetry. Actu¬ 
ally, the lepton flavor structure has been discussed in the framework of the flavor symmetry. 
Before the reactor experiments reported the non-zero value of ^ 13 , there appears a paradigm 
of ” tri-bimaximal mixing” (TBM) [321 El], which is a simple mixing pattern for leptons and 
can be easily derived from flavor symmetries. Some authors succeeded to obtain the TBM 
in the A 4 models |35]-|39]. After those successes, the non-Abelian discrete groups are center 
of attention at the flavor symmetry |in]-[l3]. The other groups were also examined to give 
the TBM IB]-[55]. The deviations from the TBM were estimated in many works [12], [56]- 
[ 74 ] . The observation of the non-vanishing 6^13 enables the detail studies of flavor models in 
the context of the sum rules of mixing angles 031, |7a-|HBl. The neutrino CP violation is 
also discussed in the context of the generalized CP symmetry [15] IM] 185] . The non-Abelian 
discrete groups with the CP symmetry have predicted the magnitude of the CP violating 
phases lEsi-iinsi. Furthermore, the neutrino mixing angles have been predicted linking with 
the quark mixing matrix by using some GUT models |106] - [TB] . Thus, the flavor models in 
the lepton sector confront the neutrino experimental data of CP violating phases. 

However, the flavor models do not always give the unique predictions since they have 
many parameters. Therefore, it is desirable to build a model with the small number of 
parameters for testing it. This attempt was proposed with Occam’s Razor [122] . where the 
inverted mass hierarchy is predicted. In this paper, we build two neutrino flavor models 
with the small number of parameters as much as possible in the framework of the indirect 
approach with the 5*4 flavor symmetry. The neutrino mass matrix is given by two complex 
parameters. The charged lepton mass matrix gives the non-vanishing reactor angle ^ 13 . The 
models lead to the inverted mass hierarchy, and predict the Dirac phase and the Majorana 
phases. Our predicted effective neutrino mass of the 0i'(3(3 is correlated with the Dirac phase. 

In section 2, we propose two simple 5*4 models, and discuss their implications. In section 3, 
we present the numerical analyses of the CP violating phases as well as the mixing angles, and 
then predict the effective neutrino mass for the Oz//?/?. The Section 4 is devoted to discussions 
and summary. In Appendices, we present the multiplication rules of the S 4 group, and the 
neutrino mass matrix in the relevant basis of neutrinos. 
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2 /S '4 flavor models 


Let us build simple lepton flavor models with the S'4 group by the indirect approach of the 
flavor symmetry. We obtain the lepton mass matrices by introducing the relevant flavons 
and assuming the alignment of their vacuum expectation values (VEVs). The advantage of 
the S'4 group includes both the doublet and the triplets as the fundamental representations. 
Therefore, the S'4 flavor symmetry was easily extended to the quark sector |107l 1108] . We 
obtain two models with the inverted neutrino mass hierarchy, which give two cases of the 
lightest neutrino mass m 3 = 0 and m 3 7 ^ 0, respectively. The particle assignments are same 
in both cases, while the vacuum alignments of flavon for the neutrino sector are different in 
each model. 

2.1 ^4 flavor model for m3 = 0 

We present a model with the S'4 flavor symmetry in the case of the vanishing lightest neutrino 
mass. The particle assignments are shown in Table [T] These assignments are similar to the 
model in Refs. jl07l 1108] except for flavon fields. Under the S'4 group, the left-handed lepton 
doublet L = (Ze, Z^, Lr) are assumed to transform as the triplet, while the right-handed 
charged lepton are assigned to the doublet Ir = hr) and the singlet tr. The right- 
handed neutrinos are also assigned to the doublet N’^ = and the singlet 

The Z4 charge is assigned relevantly to the leptons. In order to get the natural hierarchy 
between the muon and tauon masses, the Froggatt-Nielsen mechanism |123] is introduced as 
an additional flavor symmetry, where 0 denotes the Froggatt-Nielsen flavon. On the 

other hand, the Higgs doublet H is assigned to the S'4 singlet. The gauge singlet flavons are 
assigned to the triplet or triplet-prime, which have different Z4 charges as seen in Table 1. 

We can now write down the S'4 x Z4 x U{ 1 )pm invariant Lagrangian for the Yukawa 
interaction in terms of the S'4 cutoff scale A and the U{ 1 )fn cutoff scale A as follows: 

Cy = yiLlRH(j),eV{AA^) + y',llRH4>',Q^/{AK^) + y”lrRHc^y\ 

+ yolN^Hc^^/A + y'j,LN^^^Hct>'JA 

+ , ( 1 ) 

where y's are Yukawa couplings. Mi and M 2 are the Majorana masses, and H = iT 2 H*. In 
this setup, we discuss the lepton mass matrices. 
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Table 1: The assignments of leptons, Higgs and flavons. 


Let us begin with discussing the neutrino sector. In order to desirable mass matrices. 
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the relevant VEV alignments are required. We will show the potential analysis to derive the 
VEV alignments in subsection 2.3. 

We take the VEVs of the relevant flavons and VEV alignments as: 

{H)=v, (0^) = n^(l,0,0), (0'^) = <(1,1,1). (2) 


By using the multiplication rules in Appendix A, we obtain the Dirac neutrino mass matrix 
as 

/o y'nK\ ^ 

Md = 0 y'j^v'^ ^ , (3) 

Vo 0 y'^v'J 

and the right-handed Majorana neutrino mass matrix Mjy as 

/Ml 0 0 \ 

Miv = 0 Ml 0 . (4) 

Vo 0 Ms/ 


By using the seesaw mechanism, 
as 

/a + b b 

My = I ^ b 

\ b b 


the left-handed Majorana neutrino mass matrix My is given 


^ 2{yDVyvf ^ {y'pKvf 

3MiA2 ’ M 2 A 2 


(5) 


where a and b are complex parameters. In order to compare with the TBM, we move to the 
TBM base. Then, we can easily seen the family structure of the neutrino mass matrix. After 
rotating My with the TBM matrix Vtbm: 


Vtrm — 


the neutrino mass matrix turns to 


My — Vtbm 


/ ^ 

1 1 

Ve V3 





0 \ 

1 

^2 

72 ) 


-I- 36 0 < 


0 


TBM 


( 6 ) 


(7) 


This mass matrix gives us the inverted mass hierarchy with the one vanishing mass. The 
flavor mixing is deviated from the TBM only in the rotation of the (1-2) axis. 

In order to get the mass eigenvalues and the flavor mixing angles of neutrinos, we examine 
Ml My as 

/ flap ^\a\^ + \/2a*b 0\ 

Mj^My = Vtbm I \/2a6* \\al^ + Q\b\'^ + ab* + a*b 0 j ^/bm 

Vo 0 0/ 

/ ||a|2 ^|a|2 + <2|a||6|e-*^ 0\ 

= Vtbm I -|-\/2 |a|||a/-|-9|6/-|-2 |a||6| cos</? 0 j I^tbm ) (8 ) 

Vo 0 0/ 
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where a = |a|e*^“, b = and = tpa — Vb ■ Hereafter, we define |a| = a and |&| = h for 

simplicity. Then, the nentrino mass eigenvalnes are 


ml 

ml 

ml 


^ + 9b^ + 2ab cos (p — \J (a^ + 96^ + 2a&cos(p)^ — 

^ + 9b^ + 2ab cos (p + (a^ + 9b‘^ + 2abcospif — IQaP'b"^ 

0 . 


(9) 


Therefore, the nentrino mass squared differences are given as 


= ml 


ml = - 
^ 2 


+ 9b^ + 2 ab cos (p + (a^ + 9b‘^ + 2 ab cos ip)'^ — 

(a^ + 96^ + 2a6cos(p)^ — , 

and neutrino mixing matrix t/j, is 


Am^„i = ml 


U, = V, 


TBM 


COS 7] e sin rj 0^ 
cos T] 0 


sm 7] 


0 


0 


where 


tan 27] = 


2\/2a^ a? + 96^ + 6 ab cos p) 


tan-^ = 


3b sin (p 


a + 3b cos ip 


( 10 ) 


( 11 ) 


( 12 ) 


276^ — + 6ab cos ip 

Thus, the neutrino mass hierarchy is only inverted one and the neutrino mixing angles are 
determined by two neutrino mass squared differences and relative phase ip in our model. 
Next, we consider the charged lepton sector. Taking VEV and VEV alignments as0 


(0) = 0, (0^) = (na,n£2,0), (0^) = ^^(l, 0, 0), (0") = ^"(O, 0,1), 


(13) 


the charged lepton mass matrix is given as 


Mp = 



(14) 


where 


02 


= yivav- 


be = y£Vi2V 


02 


ce = y[v[v 


de = yX^j. (15) 


AA2’ AA2’ 

In order to reproduce the relevant mass hierarchy between the muon and the tauon, we take 
0/A as the Cabibbo angle 0.22 approximately. 


^We take the VEV of cj)p as {vpi,vp 2 , 0) for simplicity. As the alternative choice, we can introduce one more 
field (j)i, which charge assignment of S 4 and Z 4 is same as (j)e, and take VEV alignments as {(jii) = (f«i, 0, 0), 
i'^e) = ( 0 ,ve 2 , 0). By using these VEV alignments, we obtain the same mass matrix in Eq. (|T4)l . 
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The left-handed mixing of charged leptons is given by investigating Mi 



MiMI = 


( i(|a.P + lQp) 

\ 0 


1 

3 


{ai - y/3ci)b} 
0 


/ i(|a.P + |QP) 

\ 0 


0 

0 

\di 


liMbe\e 




2 

3 


- V^\bi\\ci\ 

\bi\^ 

0 


3 '^^hc ] 



( 16 ) 


where ai = |a£|e*‘^“^ be = \bi\e"^>‘e, a = and ip ah = Pat, - Pbe, Pbc = Pb^ - Per 

Hereafter, we dehne \ai\ = ai, \bi\ = be, and |q| = q for simplicity. Then, the charged lepton 
masses are given as 


rria g 


m,, = - 
^ 3 


2 I l2 I 2 

(i£ — 

2 I l2 I 2 I 

<2/ + 6/ + Ctf + 


j + bj + eff — bj (saj + cj + 2y/3aiCe cos 


Pe 


y {aj + bj + c]f - bj (^3aj + cj + 2V3aiCi cos 


Pe 


ml = dj , 


(17) 


where pi = pab + Pbc ■ 

The U{1)fn symmetry guarantees the mass hierarchy of the muon and the tauon. Now, 
we discuss the electron and muon masses. If we assume be S> ai, ci in Eq. (lT4ll . the electron 
and muon masses in Eq. (ITT)) are approximately written as 


ml ~ 


0a£ + + VSaiCe cos p^ , 

J - J (^26^ + - V3aic,cos:fi, 


( 18 ) 


Then, the muon mass is ■\/2/ 3bi and the electron mass is reproduced by tuning ai and q. 
On the other hand, the left-handed charged lepton mixing matrix Ui is 


cos A e sin A 0^ 
Ue = I —sin A cos A 0 

0 0 1 . 


(19) 


where 


2 y/SaiCi cos pi ai sin pab + V^ci sin pbc 

tan2A = — -—r:i-- , tan-^^ =- j= - . (20) 


bl + 


Qj ^ COS 


y/Sci 


COS ^bc 


In the charged lepton sector, there are four real parameters (a^, be, ce, de) and two phases. 
After inputting charged lepton masses, there remains two independent parameters, which are 
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the mixing angle A and the CP violating phase Those free parameters are adjusted to the 
experimental data of the lepton mixing matrix, namely Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) matrix t/pMNS |124l 1125] as 



! cos A 

—e ®®^^ sin A 

0 ^ 

i i 

f cos 77 

e ®®^ sin 77 

o\ 

t i 

f^PMNS = UjUi, = 

e®®^^ sin A 

cos A 

0 

Ptbm 

—e*^ sin 77 

cos 77 

0 

\ 

( 0 

0 

l) 

1 1 

( 0 

0 

1 / 


The mixing matrix elements are written as 
cos 7] (2 cos A + sin 


Uel = 


Vq 


in A) ^ sin rj (^— cos A + e sin A) 


x/3 


cos r] (cos A — e sin A) e sin r] cos A + sin A) 

Ue2 = - 7^ -1" 


U^i = 

U^2 = 

Url = 


\/3 V6 

cos 7] (— COS A + 2e®'^^ sin A) e®'^ sin t] (cos A + sin A) 

71 71 ’ 

cos 7] (cos A + e®®^^ sin A) e“®®^ sin r; (— cos A + 2e®®^^ sin A) 


73 

cos rj e®^ sin t] 


76 73 ’ 

e“®®^^sinA 

Ce3 = - 7^ - , = 


+ 


Ur2 = 
COS A 


COST] 


72 ’ 72 ’ 

Therefore, the mixing angles are expressed as follows: 


76 

e-iip ^ 

76 

Ut3 = —^ • 


72 


sin^ 9^■? = 




\Uel\^ + \Ue2\^ 
1 


9 + 3 cos 2A(1 + sin^ rj) + sin^ 77(1 + 8 cos ipi sin 2A) 


6(3 + cos 2A) 
cos^ 77(1 + 3 cos 2A + 8 cos ipi sin 2A) 


+ 272 sin 277| cos"0(1 + 3cos2A — cos'0£sin2A) — 3 sin 7 /) sinsin 2A 


sin^ 023 = 




\U,3? + \U^3? 


= 1 


2 - sin^ A 


sin^0i3 = \Ue3\^ = 


sin^ A 


( 21 ) 


( 22 ) 


(23) 


In order to compare our model with the TBM in the neutrino sector, it is useful to write 
down the mixing angles without the contribution from the charged lepton sector, that is ones 
in Eq. flTT]) . 


sm 


7 i 2 = ■^ (^1 + sin^ 77 + 72 sin 2r] cos -0 j , sin^ ^23 = ^ > 


sin 0 + = 0 , 


(24) 


where 9^- denote the mixing angles only in the neutrino sector. Our model is different from 
the TBM only in 0j®2. We will see how large our 9^2 is deviated from the TBM sin^ 05^2 = 1/3 
in the next section. 
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Before closing this subsection, we present a simple sum rule between 6*23 and 613 in Eq. (l2^ : 


sin^ 023 = 1 


2 COS^ 013 

By using the experimental data [126] for 3a range, the sin^ 6*23 is predicted to be 

0.487 < sin^6*23 < 0.490, 


(25) 


(26) 


which is within the 3a range for the experimental data |126] as 0.389 < sin^ 023 < 0.664. We 
show the details of our numerical analyses in section [3l 


2.2 ^4 flavor model for m3 7^ 0 

We can obtain another simple model with the non-vanishing lightest neutrino mass. The 
particle assignments are same as the model in Table [T] Therefore, the 5*4 x Z 4 invariant 
Lagrangian for the Yukawa coupling is same as in Eq. ([1]). The VEV alignments for the 
neutrino sector are different from Eq. ([2]). Taking VEV alignments as@ 

(0,) = v,{0, 1, -1) , (0^) = <(1,1,1) , (27) 

we obtain the Dirac neutrino mass matrix M^, as 

( 0 0 y'D<\ ^ 

Md = -^VdVu v'dvI j . (28) 

-^Vdv, y'nv'J ^ 

The right-handed Majorana neutrino mass matrix M^r is same in Eq. (j4]) as 

/Ml 0 0 \ 

Miv = 0 Ml 0 . (29) 

Vo 0 M 2 J 

By using the seesaw mechanism, the left-handed Majorana neutrino mass matrix M^ is 
written as 


Ml/ — \ b 2(3 b Qj b 
,b a + b 2 a + bi 


a = 


{ypVuvf 
3MiA2 ’ 


M 2 A 2 


b = 


(30) 


where a and b are complex parameters. Moving to the TBM base, the neutrino mass matrix 
is given as 

a —\/ 2 a 0 \ 

Ml/ = Vtbm I — y/ 2 a 2 a 3b 0 1 . (31) 

0 0 a) 


Mhe alignments {(j)u)lvi, = (0,1,1), (0, —1, —1), (0, —1,1) also lead to the same Mi, in Eq. (1501) . 
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In order to get the mixing angles of the left-handed Majorana neutrino, we investigate 
as 


/ 3|ap -3\/2(|a|2 +|a||6|e-*^) 0\ 

= Vtbm I -3\/2(|ap-h lall^le**^) 3(2|ap3|6p-h 4|a||6| cosy?) 0 ) , (32) 

\ 0 0 |a|7 


where a = 1016 *“^“, b = and (pa — (ft = ^ ■ Hereafter, we define |a| = a and |&| = b for 

simplicity. Then, the neutrino mass eigenvalues are 


ml 

ml 


3a^ + 3b^ + Aab cos p — \J (3a^ -|- 36^ -|- 4 a 6 cos<y 9 )^ — 

3a^ -I- 3b^ + Aab cos p + \J{3a^ -|- 36^ -|- 4a6 cos p)^ — Aa^b'^ 


ml 


= a^. 


(33) 


Therefore, the neutrino mass squared differences are 


= ml - ml = 


-a^ + 36^ - 1 - 4a6 cos p + \/{3a? + 3fe^ -|- 4a6cos<p)^ — 

I V 


= 'm .2 — = 3 \J{3a^ + 3b‘^ + Aab cos p)^ — Aa^b'^ , 

and the neutrino mixing matrix Uy is 


cos T] e sin rj 0^ 
Uy = Vtbm 1 —e^'^sinr] cost] 0 

0 0 1 , 


(34) 


(35) 


where 


tan2r7 


2 a 72 (a^ + b'^ + 2abcosp) 
-|- 36^ -|- 4a5 cos p 


b sm p 

tan'0 =-;- . 

a + b cos p 


(36) 


We can easily hnd the mass ordering with m 2 > mi > m 3 for any cos(p. Therefore, the 
neutrino mass hierarchy is only inverted one and the neutrino mixing angles are determined 
by two neutrino mass squared differences and relative phase p. 

The charged lepton mass matrix is the same as the one in the previous subsection. There¬ 
fore, the PMNS mixing matrix elements are expressed in the same form in Eq. (l22|l . We have 
the same sum rule between ^23 and 6^13 in Eq. (I25l) . 


2.3 Potential analysis and VEV alignments 

In this subsection, we present the potential analysis of flavons in the framework of the su¬ 
persymmetry with the f/(l)j? symmetry. We can generate the vacuum alignment through 
F-terms by coupling flavons to driving fields, which carry the R charge +2 under U{1)r 
symmetry. We also assign R charge -|-1 to the lepton doublets, right-handed charged leptons, 
and right-handed Majorana neutrinos. In addition to driving fields with R = 2, we introduce 
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Table 2: The assignments of the additional helds. 


a new flavon 0", which does not conple to leptons and right-handed Majorana neutrinos due 
to the U{1)fn charge. The charge assignments of the additional flavon and driving helds are 
summarized in Table |2J 

Let us write down the 5*4 x Z 4 x U{1 )fn invariant superpotential for the interaction of 
scalar helds as follow: 


wv = + 2 /?'+ Vxo'P'u'PuXo , 

where y^'^, yr^^ and y^^ are arbitrary parameters. Then, the potential V is given as 

dwy 


V = 


dXn 


— lie'll + + 0£30£3)I + \yi'o{^'ll(t^”l + 0£20£2 + | 

+ + 4>u24>i2 + 4>u3(j>l3)\ 




"2 


K3^v3/ 


+ 


^l/xo(-20m0"l + 0[.20"2 + (l>'u3(l> 


" 1 

v3) 


(37) 


(38) 


where Xq = ■Coj hO; Xo- Therefore, conditions to realize the potential minimum {V = 0) are 
given from Eq. fl38|) as 


+ 0^20a + 0a0a) “ 

+ 0a0"2 + Mk) = 0) 

1/770 + 07720^2 + </’7^30a) = O 7 

■^//xo(‘/’i720i/2 ~ 4’i/3^iy3) ~ 0; 

+ ^v2^v2 + ^v3^v3) ~ O' (39) 

One of the solutions which satishes these conditions is given as 


(0a) — (0a) ~ (0a) ~ (0a) ~ (0a) ~ (07^2) — (0i/3) — 0, 

(0U) = (0'2) = (0'.3), (0:1) = (0"2) = (0''3), (40) 

then, the VEV alignments are 


(00 = (ufl, 0), (00 = n((l, 0, 0), (0") = n"(0, 0,1), 

(00 =n,(l,0,0), (00 =n'(1,1,1), (41) 


which have been used in the subsection 2.1. We can also derive the VEV alignments of the 
subsection 2.2 in the similar discussion. 
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3 Numerical analyses 


In this section, we show the numerical analyses in our 5*4 models. We use the result of 
the global analyses of neutrino oscillation experiments [12611127[ 1128] . The 3cr range of the 
experimental data [126] for the inverted neutrino mass hierarchy are 


7.02 < 


Ami 


sol 


10“^ eW 


< 8.09, 


2.307 < 


10-3 eV^ 


< 2.590, 


0.270 < sin^ 9i2 < 0.344, 0.389 < sin^ 623 < 0.664, 0.0188 < sin^ 613 < 0.0251. (42) 


Our numerical strategy is as follows. Fixing a random value for (p of Eq. flTU]) in the region 
—71 < p <77, then, a and b are determined from the experimental data of Am^^^ and Am^^^ . 
Then, we can calculate the neutrino mixing matrix [r] and t/j) from Eqs. flTT|) and (fT2|) . The 
magnitude of A is determined in sin^ 613 of Eq. fl2^ by using the experimental data [126] for 
3cr range. Finally, choosing a random value for in —tt < ^ tt, we calculate the PMNS 
mixing matrix elements in Eq. (12ip . which are adjusted to the experimental data of Eq. fl42p . 

In our numerical studies, we predict the Dirac phase and the Majorana phases. Our 
predicted PMNS matrix in the previous section should be compared with the following con¬ 
ventional parametrization [129] including the Majorana phases: 


0 

0 A 

/ C 12 C 13 

S 12 C 13 

Si3e-®^c^\ 

/e®“ 

0 

0 \ 


0 

1 —S 12 C 23 — Ci2S23Sl3e®'^‘^'’ 

C 12 C 23 — Sl2S23Sl3e®‘^‘^®’ 

'S 23 C 13 1 

0 

gi/3 

0 

0 

gJi5x / 

\ S 12 S 23 — Ci2C23Sl3e®'^‘^®’ 

— C 12 S 23 — Sl2C23Sl3e*^'^^ 

C 23 C 13 / 

VO 

0 

1/ 


Cl2Cl3e*('^' 

(S12S23 — Ci2C23Sl3e*^'^^)e*P^““) 


Sl2Cl3e 


i(Se-p) 


i{Se—Scp)'^ 


= I (-S12C23 - Ci2S23Sl3e*'^^^)e*P'' (C12C23 - Sl2S23Sl3e*'^^^)e*P'' 


Sise 
S23Cl3e 


( —C12S23 — Sl2C23Sl3e*‘^'^^)e*P^ C23Ci3e 


(43) 


where Cij and Sij denote cos 9ij and sin6*jj, respectively. The phases 5e, and 5r could be 
absorbed in the left-handed charged lepton fields, 5cp is the Dirac phase, and a, (3 are the 
Majorana phases. 

We can calculate the Dirac phase 5cp as follows. The Jarlskog invariant |13U| . which is 
the measure describing the size of the CP violation, is given as 


Jcp = Im [c,u,jU;,vp 


sin 2A 
48 


3\/2 cos 'tpi sin '0 sin 2 r] — sin I \pi cos sin 27] + A cos 2 ri 


(44) 


where x, y = e, ]i, t and i, j = 1-3. Then, the CP violating Dirac phase is written in terms 
of the lepton mixing angles and Jcp as 


Jcp 


sinhcp = — 

■S23C23'Sl2Cl2Sl3Cf3 

In our model, the one of the PMNS mixing matrix elements Uri is written as 


(45) 


IVi 


cos 7 ] sin t] 

ys 


S12S23 - Ci2C23Si3e*^°P 


(46) 


10 


























then, the CP violating Dirac phase is also written as 

+ c? 2 C^ 3 S ?3 - I cos 2 T] - ^ sin^ r; - 1 sin 2t] cos V' 
cos5cp =- o, , ^ ^ , - • (47) 

^■Sl2S23Ci2C23Si3 

Thus, 6 cp is determined up to the quadrant. 

Next, let us calculate the Majorana phases. In our model, and Urs have no phases as 
seen in Eq. (|22il . so we find = 5,- = 0 in Eq. (03]). For convenience, expressing the phases 
of the mixing matrix of Eq. (|2^ as 


Del = \Uel\e^^\ 

De2= |De2|D'^ De3 = | De31 e-'^^ 

(48) 

we derive the following relations by comparing with Eq. (05]): 


61 = Se — a , 

62 — 6 e — (3 , — ^CP — 6 e ■ 

(49) 

Eliminating 6 e in these relations, we 

obtain 


><7 

1 

1 

tt, 

0 

(3 = 6 cp — — 62 -, a — (3 = 62 — 5i, 

(50) 


from which we can calculate the Majorana phases numerically. 

3.1 <S '4 flavor model for m 3 = 0 

In this subsection, we show the numerical re¬ 
sults in the case of the lightest neutrino mass 
m 3 = 0. In order to reproduce the relevant 
mass squared differences, and Amg^j 

for the inverted mass hierarchy, a ~ 36 and 
ip ~ ± 7 r should be satisfied. Then, the -0 of 
Eq. (IT^ is restricted in 45° < It/^l < 90°. On 
the other hand, sin rj in Eq. (1121) is restricted in 
—0.70 < sin ?7 < —0.38. In FigH] we show the 
allowed region on -^-sin^ 9’(2 plane, where 612 is 
the mixing angle of the neutrino sector with¬ 
out the charged lepton contribution as seen in 
Eq. fl2T|) . As seen in this figure, our model is 
possibly deviated from the TBM considerably. 

Then, the mixing angle and the phase of the charged lepton sector become important to 
adjust the observed PMNS mixing angles. 

Now we can calculate the Dirac phase and the Majorana phases. In Figs. [2] and [S] we show 
the Dirac phase (5cp versus the Majorana phases a and /3, respectively. The both positive 
and negative signs are allowed for a and (3. The allowed regions of the Majorana phases are 
110° < |a| < 150° and 145° < |/S| < 170°. On the other hand, the Dirac phase is still allowed 
in the all region as —vr < 6 cp < tt. 

Since the difference between the Majorana phases a and 13 contributes to the 0z//9/3, we 
show the Dirac phase versus (a — (3) in Fig. 0] We obtain 10° |a — /5| ^ 50°. We also show 



Figure 1: The allowed region on -^-sin^ ^ 5^2 
plane. The red line corresponds to the TBM. 
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Figure 2: The predicted region on a-5cp 
plane. 


Figure 3: The predicted region on pScp 
plane. 



150 
100 
50 
0 

-50 
-100 
-150 

-150 -100 -50 0 50 100 150 

o-n 



Figure 4: The predicted region on (a—/?)- 
Scp plane. 


Figure 5: The predicted region on a- 
f3 plane. The Majorana phases are re¬ 
stricted in the Conner on this plane. 


the allowed region of the Majorana phases on the a-f3 plane in Fig. O The Majorana phases 
a and f3 are restricted in the two regions on this plane. 

Based on the numerical results of the phases, we can estimate the effective neutrino mass 
for the 0u(3(3. The effective neutrino mass is written as 

3 

mee = ^miUl. (51) 

i=l 

We show the effective neutrino mass of the Oi^/3/3 versus the difference of the Majorana phases 
and the Dirac phase for in Figs. Elandl?! respectively. For the hxed effective neutrino mass 
for the 0z//5/3, there are two fold degeneracy and four fold degeneracy in the Majorana phase 
(a — (3) and the Dirac phase Scp, respectively. Therefore, the degeneracy will be solved if 
both 6 cp and rriee are observed. At present, we predict 32 meV < \mee\ ^ 49 meV, which is 
close to the expected reaches of the coming experiments of the 01/(3(3 m- 
Finally, we show the sum of the neutrino masses is predicted as 

0.0952 eV < < 0.101 eV. (52) 
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-40 -20 0 20 

o'-js n 


40 


Figure 6 : The predicted effective neutrino 
mass for the versus {a — (3). 


-150 -100 -50 0 50 100 150 

ScF [d 

Figure 7: The predicted effective neutrino 
mass for the Qi'(3(3 versus 5cp- 


This total neutrino mass is within the reaches of the future cosmological and astrophysical 
measurements, such as the CMB spectrum, the galaxy distributions, and the red-shifted 21cm 
line in the future Ea. 


3.2 S/i flavor model for ms ^ 0 

In this subsection, we show the numerical re¬ 
sults in the case of the non-vanishing light¬ 
est neutrino mass m 3 7 ^ 0. In order to re¬ 
produce the relevant mass squared differences, 
Am^^rn inverted mass hierar¬ 

chy, a ~ 6 and ip ~ ivr are satisfied. Then, the 
if} of Eq. fl36ll is restricted in 50° < |'0| ^ 90° 
numerically. On the other hand, sinr^ in Eq. 
fl36|) is restricted in —0.70 < sinr^ < —0.54. 

In Fig. IHl we show the allowed region on '0- 
sin^ 6'(2 plane, where is the mixing angle of 
the neutrino sector without the charged lepton 
contribution as seen in Eq. fl2T]) . This result is 
very similar to the case of m 3 = 0 of Fig. [T]in 



Figure 8: The allowed region on -^-sin^ 05^2 
plane. The red line corresponds to the TBM. 

the previous subsection. 
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Figure 9: The predicted region on a-5cp 
plane. 




an 


Figure 10: The predicted region on 13-6cp 
plane. 
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Figure 11: The predicted region on {a 
(3)-6cp plane. 


Figure 12: The predicted region on a- 
f3 plane. The Majorana phases are re¬ 
stricted in the Conner of this plane. 


At hrst, in Figs. [9] and [TOl we show the Dirac phase (5cp versus the Majorana phases a 
and /9, respectively. The both positive and negative signs are allowed for a and f3. The 
allowed regions of the Majorana phases are 105° < |a| ;< 135° and 145° < |/3| < 170°. The 
Dirac phase is still allowed in the all region as —vr < 6 cp ^ tt. We also show the Dirac phase 
6 cp versus {a — (3) in Fig. [TT] We obtain 10° < \a — I3\ < 60°. In Fig. [T21 we show the 
allowed region of the Majorana phases a and (3. The Majorana phases are restricted in the 
two regions on this plane as well as the case of the previous subsection. 

Next, we also show the effective neutrino mass ruee for the versus the Dirac phase 

in Fig. [131 For the hxed effective neutrino mass for the Quf3f3, there are four fold degeneracy 
in the Dirac phase 6 cp as well as the case of the previous subsection. This degeneracy will 
be solved if both 6 cp and niee are observed. The predicted magnitude of niee is 34 meV < 
\^ee\ ^ 59 meV, which is close to the expected reaches of the coming experiments of the 
Qul3f3 [m] . 

Finally, we show the allowed region of neutrino masses m 2 and m 3 in Fig. [TTl The neutrino 
mass m 2 is 58.7 meV < m 2 < 62.2 meV. The lightest neutrino mass m 3 is 31.0 meV < m 3 < 
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35.8 meV. The sum of the neutrino masses is predicted as 


0.150 eV < Y^rui < 0.160 eV. (53) 

i 

This value is also within the reaches of the future cosmological and astrophysical measure¬ 
ments. 




0.055 

.. -. •••■ ..•■• 

> 0.050 

.ft . • 

* . •'•i 4*^' . ** 


.* • i. •f’ • * • • 

Ji 0.045 

¥ -r- 


■ . ■’S 

0.040 

'■C' 't'“. ..l.r 



0.035 



-150 -100 -50 0 50 100 150 


Scp [°] 

Figure 13: The prediction on Scp-\mee\ 
plane. The allowed region of effective neu- 


trino mass for the 0i^/3/3 is 34 meV < 
Irrteel < 59 meV. 


0.0620 
0.0615 
^ 0.0610 
^ 0.0605 

S 0.0600 

0.0595 
0.0590 

0.0340 0.0345 0.0350 0.0355 

m 3 [eV] 

Figure 14: The allowed region on m^-m 2 
plane. The lightest neutrino mass m 3 is 
proportional to m 2 . 



4 Discussions and Summary 

We have presented the minimal 5*4 flavor models of leptons. By using the alignments of the 
VEV’s of the flavons obtained in subsection 2.3, the models lead to the two different neutrino 
mass spectra with the inverted mass hierarchy. The neutrino mass matrix is given by two 
complex parameters in Eqs.([5]) and (l30D . The hrst case has one vanishing neutrino mass 
m 3 = 0, and the second case has three non-vanishing neutrino masses. The charged leptons 
are not flavor diagonal. There remains the hrst and second family mixing, which give the 
non-vanishing reactor angle ^ 13 . We have moved to the TBM basis of neutrinos in order to 
clarify the deviation from the TBM. 

It is not necessary to go to the TBM basis to calculate the mixing angles, but this is 
merely a convenient basis to do so. We present alternative basis, the bi-maximal base or the 
p-r symmetric one, where it is also convenient to study the mass matrices in appendices B 
and C, respectively. 

Inputting the experimental data of two neutrino mass squared differences and the mixing 
angles 6 ^ 12 , ^23 and ^ 13 , we can predict the Dirac phase and the Majorana phases. Furthermore, 
we predict the magnitude of the effective neutrino mass of the 0i^/3/3, which is correlated with 
the Dirac phases. For the hxed effective neutrino mass for the 0z//9/9, there are two fold 
degeneracy and four fold degeneracy in the Majorana phase difference {a — (3) and the Dirac 
phase 6 cp, respectively. Therefore, if both 6 cp and mge are observed, the Majorana phases 
are determined hnally. At present, the predicted magnitudes of rriee are in the regions as 
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32 meV < |mee| < 49 meV and 34 meV < \mee\ 59 meV for the cases of m 3 = 0 
and m 3 7 ^ 0, respectively. These values are close to the expected reaches of the coming 
experiments of the 

The sum of three neutrino masses is also predicted in both models as 0.0952 eV < ^ m* < 
0.101 eV and 0.150 eV < ^m^ < 0.160 eV, respectively. Those total masses are also within 
the reaches of the future cosmological and astrophysical measurements. 
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Appendix 

A Multiplication rules of the S4 group 

We show the multiplication rules of S 4 : 




( ai 
02 
03 


3 


0 



= {aibi + 0262)1 © (-0162 + 0261)1' © ^ 

/ 0261 \ / ai 6 i \ 

= I -i(\/ 3 ai 62 + 0262) I ©I |(\/ 3 a 262 - 0162) 1 

\ |(\/3ai63 - 0263 ) / 3 \-|(\/3a263 + 0163 )/ 3, 

( 0161 \ / 0261 \ 

i(\/ 3 a 262 - 0162) 1 © I-|(\/ 3 oi 62 + 0262) 1 

-i(\/3a263 + 0163 )/ 3 \ |(\/3ai63 - 0263 ) / 3, 


/ 7 7 7 \ -^(0262 — 0363) \ 

= (“li’i + <7262 + < 7363)1 e + aA))^ 

/ 0263 + 0362 ^ / 0362 - 0263 

© 0163 + 0361 © 0163 - 0361 

\ai 62 + 0261 / 3 \a 26 l - 0162 

/ 7 7 7 \ -^(0262 — O 363 ) \ 

= (a,61 + ajfti + 0363)1 e + a,b2 + 0363))^ 

( 0263 + 0362 ^ /0362 - 0263 

O 163 + 0361 j © I 0163 — 0361 

0162 + (l2bl) 3 \O261 — O162 




— (O 161 + 0262 + 0363 ) 1 ' © 



(2ai6i — 0262 — 0363) 
^(0262 - 0363) 


2 


/0362 - 0263 
© j 0163 - 0361 
\a 261 - 0162 


3 


/ 0263 + 0362 ^ 

I 0163 + 0361 
\a 162 + 0261 / 3 , 


More details are shown in the review [HI 02]. 


(54) 
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B Neutrino mass matrix in the bimaximal base 


B.l m 3 = 0 case 

The left-handed Majorana neutrino mass matrix is 



After rotating My with the bimaximal mixing matrix Vbm as 

0 \ 


l^BM — 


/— — 

v/2 

~2 2 

_ 1 1 

2 2 


1 

V2 / 


the left-handed Majorana neutrino mass matrix is rewritten as 


My — 1/bm 


'l{a+{3-2V2)b) l{a-b) 


0^ 


( 55 ) 


(56) 


0 


b) I (a -|- (3 -|- 2^/2)b) 0 j I^bm 

0 0, 


(57) 


B.2 m-i ^ 0 case 

The left-handed Majorana neutrino mass matrix is 

(h h b 
Mjy “ j ^ T ^ Cl b 
\b Cl b 2(2 bj 


(58) 


After rotating My with the bimaximal mixing matrix I^bm, the left-handed Majorana neutrino 
mass matrix is rewritten as 


'I (3a + (3 - 2V2)b) 
My = Vbm | —|(3a -|- b) 

0 


— 4 (3a -|- 6) 


0^ 


I (3a -|- (3 -|- 2\/2)5) 0 j Idgjyj . 
0 


(59) 
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C Neutrino mass matrix in the /i-r maximal base 

C.l m3 = 0 case 

The left-handed Majorana neutrino mass matrix is 

+ h b 


M, = 


b b b 
b b bj 


After rotating My with the /i-r maximal mixing matrix V 23 as 


A 0 


0 


^ 23 - (0 ^ ^ 

v/2 

the left-handed Majorana neutrino mass matrix is rewritten as 

/a + b \/2b 0 \ 

My = 1^23 \ V2b 2b 0 Vi . 
\ 0 0 0 / 

C.2 m3 7^ 0 case 

The left-handed Majorana neutrino mass matrix My is 


My — \ b 2(3 b CL b 
yb CL b 2(3 T 

After rotating My with the /i-r maximal mixing matrix V 23 , the left-handed Majorana 
trino mass matrix is rewritten as 


V2b 0^ 


rT 


My = V23 {V2b 3a + 26 0 | . 

0 0 a. 


( 60 ) 


(61) 


(62) 


(63) 

neu- 

(64) 
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